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Abstract 

We present an ab initio theoretical formalism for the static paramagnetic spin sus- 
ceptibility of metals at finite temperatures. Since relativistic effects, e.g. spin-orbit 
coupling, are included we can identify the anisotropy or easy axes of the spin fluctua- 
tions. Our calculations find hep-iron to be unstable to in ab-plane, incommensurate 
anti-ferromagnetic (AFM) modes (linked to nested Fermi surface) below T/v =69K 
for the lowest pressures under which it is stable. T\r swiftly drops to zero as the 
pressure is increased. The calculated susceptibility of yttrium is consistent with the 
helical, incommensurate AFM order found in many rare-earth-dilute yttrium alloys. 
Lastly, in line with experimental data, we find the easy axes of the incommensu- 
rate AFM and ferromagnetic spin fluctuations of the normal state of the triplet 
superconductor Sr2Ru04 to be perpendicular and parallel with the crystal c-axis 
respectively. 
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At low temperatures under sufficient pressure iron swaps its body-centered cubic crystal struc- 
ture in favor of a hexagonal close packed arrangement with the loss of its ferromagnetic long range 
order. It also becomes superconducting [1] in a narrow range of pressures between 15 to 30 GPa. 
This recent observation of superconductivity has added the phase of a common element to the grow- 
ing list of magnetic or nearly magnetic materials exhibiting a superconductivity in which magnetic 
fluctuations seem to play a key role [2]. Prior to this finding the main interest in the magnetic prop- 
erties of the high pressure hexagonal phase of iron was motivated by their effect upon its structural 
and mechanical properties. As part of the study of the phase diagram of iron this has a wide spread 
relevance for subjects as disparate as, for example, the properties of steel and the structure of the 
Earth's inner core [3]. 

The importance of magnetism upon the phase diagram is well-known. Ab-initio density func- 
tional theory (DFT) calculations have demonstrated that the formation of magnetic moments sta- 
bilises the b.c.c. phase of iron at low temperatures and normal pressures. Although stable only 
above the Curie temperature, face-centered cubic iron has been found to have an incommensurate 
spin density wave (SDW) magnetic state at low temperature following measurements on fcc-Fe pre- 
cipitates on copper [4] and ab-initio calculations [5] which also show the magnetic interactions to 
vary sharply function of volume. 

The magnetic attributes of hexagonal iron have not been so well identified. To date DFT calcu- 
lations of the energies of two possible anti-ferromagnetic (AFM) ground states are lower than the 
non-magnetic one for pressures up to 60 GPa and calculated elastic constants for such magnetic 
states provide a better interpretation of experimental data than those from non-magnetic states [6]. 
On the other hand, in-situ Mossbauer data from hep iron [7] show no evidence of long range mag- 
netism although they do allow for it to possess strongly enhanced paramagnetism which diminishes 
with increasing pressure. In this letter we study the spin fluctuations (SF) in this material via 
calculations of the paramagnetic spin susceptibility as a function of both temperature and pres- 
sure. We find the predominant modes to be incommensurate with the hexagonal lattice and set 
by nesting features of the Fermi surface. We find their 'easy axis' to be in the ab-plane. These 
fluctuations become unstable at low temperatures and at the lower pressures where hcp-Fe is just 
stable, indicating an incommensurate SDW magnetic ground state. 

To carry out this study we describe a theoretical formalism to calculate a relativistic param- 
agnetic static spin susceptibility for metals at finite temperature. We highlight how the natural 
inclusion of spin-orbit coupling allows us to study the anisotropy of enhanced spin fluctuations 
in nearly magnetic systems and to identify their easy- axes. The paramagnetic susceptibility of a 
metal can possess much structure in wave-vector (q) space. Indeed a peak at one temperature 
can evolve into a divergence at a lower one to signify the system's transition into a magnetically 
ordered state characterised by a static magnetisation wave with a wave-vector corresponding to 
that of the peak. By seeking out the peaks of the paramagnetic susceptibility by examining a wide 
range of wave-vectors we can identify the dominant spin fluctuations and probe for many differ- 
ent possible magnetic ground states. Moreover the easy axis of the enhanced spin fluctuations in 
the paramagnetic phase will indicate the direction along which magnetisation grows in the ordered 
phase. 

Following the application to hcp-Fe, we calculate the q-dependence and easy axis of spin- 
fluctuations in hep metals yttrium and scandium. Our results are consistent with the helical 
anti-ferromagnetic structure that occurs when Y is doped with magnetic rare earth impurities 
such as Gd. Finally we investigate the anisotropy of the spin fluctuations in the normal state of the 
perovskite spin-triplet superconductor Sr 2 Ru04. This aspect has been proposed as an ingredient 
for the theoretical understanding of the pairing mechanism in this extensively studied material [8]. 



Our principal findings are in good agreement with NMR data [9]. 

We access the paramagnetic spin susceptibility by considering a paramagnetic metal subjected 
to a small, external, inhomogeneous magnetic field, 5h ext '(r), which induces a magnetisation 5m(r). 
We use DFT to derive an expression via a variational linear response approach [10,11]. Although 
there are a number of studies of this type [12], relativistic effects are typically omitted. We address 
this issue in this letter. Recently we have developed [13] a scheme for calculating the wave-vector, 
frequency and temperature-dependent dynamic spin susceptibility and have applied it to Pd, Cr 
and Cr-alloys obtaining good agreement with experimental data. Here we describe the new aspects 
that emerge when relativistic effects are incorporated in its static limit. 

We start with the relativistic version of density functional theory [14] where a Gordon decom- 
position is applied to the current density and the spin-only part of the current retained. The Kohn- 
Sham-Dirac equation involving effective one-electron fields for a paramagnetic system is solved by a 
one-electron Green's function G (t, r'; e). If a small external field 5h ext is applied along a direction 
n with respect to the crystal axes, a small magnetisation, 5m(r), and effective magnetic field are 
set up. The latter is <5b e ^[p(r), m(r)] = 5h ext (r) + I xc (r)5m(r) where we have assumed that Sb ext 
couples only to the spin part of the current and I xc (r) is the functional derivative of the effective 
exchange and correlation magnetic field (within the local density approximation e.g. [15] LDA) 
with respect to the induced magnetisation density. The response of the system to the external 
magnetic field is expressed in terms of the Green's function obtained from a Dyson equation in 
terms of unperturbed Green's function G (r,r';e) of the paramagnetic system and 5b e f f . For a 
general crystal lattice with N s atoms located at positions ( I — 1,..,N S ) in each unit cell, a lattice 
Fourier transform can be carried out and an expression for the full static spin susceptibility found. 

X n (x,, xj, , q) = x^l, 4 , q) + £ / Xo(*i, x£, , q) 4c(xf„) jftf,,, xj, , q) dx£, (1) 
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with the non-interacting susceptibility of the static unperturbed system given by 
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The integral is over the Brillouin zone with wave vectors k, q and k + q within the Brillouin zone 
of volume vbz- The x; are measured relative to the positions of atoms centred on a;. The fermionic 
Matsubara frequencies u n are (2n + l)7r/c B T and // is the chemical potential. The Green's function 
for the unperturbed, paramagnetic system containing the band structure effects is obtained via 
relativistic multiple scattering (Korringa-Kohn-Rostoker, KKR) theory [16]. The susceptibility has 
a dependence on the direction of the magnetic field, n, which is lost when relativistic, spin-orbit 
coupling effects are omitted. If an external magnetic field is applied along the (0, 0, 1) direction, 
i.e. the c-axis in h.c.p. and tetragonal systems, we have \ z whereas \ x is given for the field in 
the ab-plane, n = (1,0,0). We obtain an anisotropy as the difference in the susceptibility when 
an external magnetic field is applied in two distinct directions with respect to the crystal axes i.e. 
(x x — X z )- This approach presented here is applicable to ordered compounds and elemental metals 
and can be modified to study disordered alloys [13]. 

Equation (1) is solved using a direct method of matrix inversion. The full Fourier transform 
is then generated x A (q, q) = (1/V) £j e iq - (a '~ a i' ) / dx ; / dx! v e^ 1 ^ X h (^i, 4 , q) where V is 
the volume of the unit cell. Some aspects of the numerical methods used to evaluate (2) and (1) 
can be found in [13,17] and further details are given elsewhere [18]. We have used atomic sphere 
approximation (ASA), effective one-electron potentials and charge densities in the calculations for 
hcp-Fe, Y and Sc. For Sr 2 Ru04 we have used muffin-tin (MT) (1=0) components of the potentials 



from a full potential calculation [19]. In all cases the details of the electronic structures including 
the Fermi surfaces compare well with those from full potential calculations. 

Figure 1 shows some results of our susceptibility calculations for hcp-Fe for four volumes in the 
range (128oq < V < 143(Iq) where ao is the Bohr radius. The largest volume we consider is 142.23 
which corresponds to a pressure of P~16GPa [20] and the smallest volume is 128.98 Oq or P~45GPa. 
For all the volumes the axial ratio of the lattice constants is taken to be c/a = 1.6 which is very 
close to the experimentally observed value over a wide pressure range [20]. Calculations of x(q) 
were performed for many wave-vectors in the Brillouin zone in a thorough search for dominant spin 
fluctuations and potential magnetically ordered states. The important ones are shown in Fig.l. The 
temperature range for such calculations was 50K < T < 300K. According to the full potential DFT 
calculations of Steinle- Neumann et al [6], two AFM configurations (termed AFM-I and AFM-II) are 
more stable than non-magnetic or FM configurations for pressures up to ~ 60GPa at T = OK. AFM- 
I has magnetisation alternating in orientation in ab-planes stacked along the c-axis while AFM-II 
has magnetisation with opposite orientation on each layer aligned with the c-axis and perpendicular 
to the x-axis. If AFM-I were the ground state magnetic configuration then for temperatures above 
the Neel temperature, T/v we should see a peak in the susceptibility at the special point wave- vector 
°La — (0, 0,a/2c), (in units of 1-nja). Similarly if AFM-II was the ground state configuration then 
we expect to see the peak in x(q) a ^ the special point qM = (l/\/3, 0, 0). For all four volumes, 
however, we find the maximal peak in x(q) to be at an incommensurate wave- vector lying in the 
basal plane of q iriC . = (0.56, 0.22, 0). We can trace this wave-vector q iriC . to a nesting of the metal's 
Fermi surface as shown in Fig.2. This is a cross-section of the Fermi surface which is generated from 
4 bands straddling the Fermi energy and is comprised of four sheets. A cross-section through the 
basal plane reveals the Fermi surface to be dominated by two hexagonal- like shapes centred on T. 
The incommensurate qj nc .=(0.56, 0.22, 0) nests two pieces of Fermi surface as shown in Fig.2 which 
is a cross-section where the x-axis is along the direction of qj nc . and the y-axis is along the c-axis. 

These results suggest that for the larger volumes the ground state of hcp-Fe is an incommensurate 
AFM. As the pressure is increased the Neel temperature decreases to zero. Figure 1 shows that 
at 100K enhanced AFM-SF exist which decrease in strength as volume decreases with increasing 
pressure. Further calculations of the full spin susceptibility, eq.(l), reveal that an instability arises 
at the nesting vector q inc . above 50K at T N = 69K for the volume V = 142.23 (a|j) (P~16GPa) where 
hcp-Fe is stable. T N drops below 50K (the lowest temperature we can consider) for the other smaller 
volumes. We therefore infer that hcp-Fe has an incommensurate SDW ground state in a small 
pressure range starting at the onset of the hep-phase ~15 GPa. Our calculations of the anisotropy 
reveal xliOinc) > xli^inc.) for all 4 volumes and all temperatures in the range (50K<T<300K). 
Thus hcp-Fe exhibits incommensurate AFM-SF in the ab-plane. If the superconductivity of hcp-Fe 
is unconventional as has been suggested [2] then the in-plane incommensurate AFM-SF must play 
an important role in the pairing mechanism. 

Fig.l also shows how x _1 (q) varies with volume for other salient wave-vectors at T = 100K. We 
see that x(<knc.) > x{<Ik) > x(qiw) > x(Qa) > x(q - (0,0,0)) for all 4 volumes. Apparently hcp- 
Fe is furthest away from forming a ferromagnetically ordered state as shown by the relatively high 
value of x _1 (qo)- This is followed by AFM-I (q^) and then AFM-II (qM)- In agreement with [6] we 
find that the AFM-II is a more stable configuration than AFM-I but we find a third special-point 
antiferromagnetic structure, AFM-III, characterised by q^ = (1/^/3, 1/3, 0), to be more stable than 
either of these. This structure would have the magnetisation direction alternating in layers stacked 
along q K . Although q K signifies a more stable AFM configuration than either AFM-I and AFM-II, 
it is at the incommensurate vector (qj nc .) where we see the maximal peak in x(q) which leads to 
our prediction that hcp-Fe has an incommensurate SDW ground state. 



Having established a way of determining the easy axis for spin fluctuations with our relativistic 
formalism, we can use it to deduce how magnetic impurities are likely to orientate in paramagnetic 
host metals. Excellent case studies are provided by yttrium alloys with low concentrations of rare 
earths. When crystal field effects of the rare earths do not dominate, the enhanced susceptibility 
of Y determines the magnetic ordering. Here we describe our calculations of this for Y and also its 
lighter 3d counterpart Sc. Both metals have hep crystal structures with very similar electronic band 
structures and Fermi surfaces. We have used experimental lattice constants a = 6.89, c = 10.83 for 

Y and a = 6.24, c = 9.91 for Sc, in atomic units, ao- Our calculations of the static spin susceptibility 
defined by (1) show a peak at the wave-vector q inc = (0, 0, 0.18)27r/a = (0, 0, 0.577r/c) for both Y 
and Sc. This vector q iriC once again corresponds to a nesting effect from a webbing feature of 
flat parallel sheets [21] as confirmed by Fermi surface calculations and is in good agreement with 
experiments [21,22]. Fig. 3 shows the temperature dependence of the anisotropy of Y and Sc at 
the nesting vector qj nc . and also at a small vector q G ~ (0, 0, 0) (indicating ferromagnetic spin 
fluctuations (FM-SF)). Evidently the anisotropy is an order of magnitude larger in Y that in Sc 
- e.g. at T = 50 K the anisotropy (Xoi^nc.) — Xo(<W.)) for Y is 0.0411 (fi 2 B /Ry.) while for Sc it 
is 0.00635 (fi B /Ry.). This difference is caused by the spin-orbit coupling being more pronounced 
in the heavier 4d metal Y than in the 3d Sc. At qj nc ., the AFM-SF have an easy direction in 
the ab-plane, xS(Qmc) > xliQinc), whereas that of the FM-SF is along the c-axis, Xo > Xo- The 
anisotropy is also larger at qj nc . than for small wave-vectors. These calculations show that both 

Y and Sc exhibit AFM-SF characterised by the incommensurate wave- vector q inc = (0, 0, 0.57tt/c) 
and are aligned in the basal plane. They explain why when Y is perturbed by the addition of rare 
earth magnetic impurities (e.g. Gd), it typically responds by ordering the impurity moments into 
a helical AFM state in which the moments align in the basal plane and rotate their orientations in 
successive planes around the c-axis [23]. 

We conclude with our study of the anisotropy of the spin-fluctuations of the normal state of 
Sr 2 Ru0 4 . We have used lattice constants of a = 7.30, c = 24.06 in a for the tetragonal unit cell. Our 
calculations of the susceptibility and Fermi surface for Sr 2 Ru04 show an incommensurate nesting 
vector of qj nc .=(0.35, 0.35, 0) similar to that found in previous calculations [24] and experiment [25]. 
We have compared the anisotropy at the incommensurate wave-vector qj„ c with that at small wave- 
vectors describing ferromagnetic SF. For the former we find (xl(<knc.) > Xoi^nc.)) implying that the 
incommensurate AFM-SF should have an easy direction along the c-axis. Conversely for small q we 
find (Xo > Xo) an d therefore the FM-SF to be in the ab-plane. These findings are in good agreement 
with NMR experiments [9]. At T = 50 K, the anisotropy (Xo(*W;.) _ Xoi^nc.)) is 0.847 (fi B /Ry.) 
which is an order of magnitude greater than what we find for Y owing to the large tetragonal 
distortion (c/a = 3.3) of the unit cell. Finally a calculation of {xZi^inc) I xli^hnc)) at T=50K yields 
0.98. According to a scenario based on a simple model by Kuwabara et al [8], the spin-triplet 
superconducting state of Sr 2 Ru04 can be stabilised by incommensurate AFM-SF with the easy axis 
for the fluctuations along the c-axis. Our value of the anisotropy of the SF, (Xoi^inc.) / Xli^inc.)) , of 
0.98, however, would also imply that the system must be very close to ordering magnetically. 
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FIG. 1. The inverse susceptibility of hcp-Fe at T = 100K for various wave-vectors (q ~ 0, the nesting 
vector q^nc. = 27r/a(0.56, 0.22, 0) and at the special points A, M and K) for a number of unit cell volumes. 




FIG. 2. Cross-section of the Fermi surface of hcp-Fe at V=137.72 a|j 
Qinc. = 27r/a(0.56, 0.22, 0) and the c-axis (0,0, 1) respectively. 
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FIG. 3. The temperature dependence of the anisotropy Xo(<l) ~~ Xo(*l) °f Y and Sc at both the nesting 
vector qinc. = (0, 0, 0.577r/c) and q ~ (0,0,0). 



